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QUANTUM  MECHANICAL  CALCULATION  OF  LIGHT  METAL  DIATOMIC 
HYDRIDES,  OXIDES,  AND  HALIDES 

ABSTRACT 

An  abbreviated  discussion  is  given  of  certain  quantum  mechanical 
calculations  completed,  in  course,  or  beginning  for  small  molecules. 
The  development  of  theory,  computer  programs,  and  extensive  calcula¬ 
tions  on  atoms  is  outlined.  This  includes  polarizability  studies, 
relativistic  calculations  and  extensive  configuration  inter  tion  for 
a  few  atomic  systems. 

Results  for  diatomic  molecules  in  minimal  LCAO-MO-SCF  approxima¬ 
tion  is  reviewed,  and  the  gradual  extension  of  this  effort  first  to 
double-zeta  LCA.C-MO-SCF  results  and  presently  to  true  Hartree-Fock 
accuracy  by  the  matrix  method  is  discussed.  The  range  of  diatomic 
molecule-systems  treated,  the  related  molecular  properties  calculated, 
and  physical  or  chemical  interpretation  effort)?  ore  outlined.  A  gen¬ 
eralized  Hartree-Fock  method  has  been  tried  out  on  Lig  and  H^  and 
promises  considerable  Improvement.  A  resume  Is  also  given  of  certain 
polyatomic  calculations  carried  out.  This  latter  effort  includes 
limited  (or  minimal)  basis  set  calculations  for  certain  linear  tri- 
atomic  molecules,  methane,  and  one-center  expansion  calculations  for 
A^  hydrides,  all  within  the  self-consistent-field  procedure.  Calcu¬ 
lations  on  the  hydrogen  molecule  using  a  lengthy  James-Coolidge  ex¬ 
pansion,  have  been  carried  to  within  a  few  kaysers  of  experiment  for 
the  ground  and  certain  few  excited  states. 
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A.  THEORETICAL  DEVELOPMENTS  AND  CALCULATIONS  ON  ATOMS  * 

The  theoretical  formulation  of  the  matrix  solution  of  the 
Hartree-Fock  equations  for  open-shell  configurations  was  studied 
extensively  in  this  period.  The  generality  of  the  allowed  open- 
shell  cases  was  progressively  extended  or  simplified  and  machine 
programs  (in  fact  a  sequence  of  such  programs)  were  constructed  to 
carry  out  such  calculations.  The  theory,  numerical  methods,  and 
programming  were  unusually  closely  correlated  as  a  unit.  An  ex¬ 
tensive  table  of  open-shell  vector-coupling  coefficients  was  also 
obtained  for  cases  involving  several  open-shells  of  different 
symmetries.  Using  these  programs,  extensive  calculations  of  wave 
functions  and  energies  to  Hartree-Fock  accuracy  were  obtained  for 
first  and  second  row  atoms  (in  collaboration  with  others) .  Exten¬ 
sive  calculations  were  also  carried  out  for  ls22sm2pn  states  of  N, 

N  ,  0,  and  0  ;  ground  and  several  excited  states  of  A£,  Cu  ,  and  Cu 
ground  states  of  Mo  and  Mo+;  and  the  ground  states  of  K+,  K,  and  Ca 
A  careful  and  extensive  study  of  the  K-shell  and  X-ray  transition 
probabilities  of  some  Neon-like  and  Argon-like  ions  was  carried  out 
using  calculated  SCF  wave  functions.  Still  in  progress  are  calcula 
tions  on  transition  probabilities  and  wave  functions  for  certain 
relevant  states.  Complementary  to  tnis  effort,  and  continuing, 
is  the  programme  to  numerically  sol\e  the  Hartree-Fock  equations 
for  atoms .  This  latter  effort  has  led  to  the  development  of  newer 
and  superior  numerical  methods  for  this  old  problem. 
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The  fully-coupled  self-consistent-field  equations  for  a  closed 
shell  atom  in  an  electric  field  have  also  been  formulated,  coded, 
and  solved.  The  dipole  polarizabilities  of  He-like,  Be-like,  Ne-like, 
and  Mg-like  ions  have  thus  been  calculated  in  good  agreement  with 
experiment  (where  known,  except  for  H~) . 

A  relativistic  self-consistent-field  theory  has  been  developed 
for  closed  shell  atoms.  In  this  theory  the  many-electron  Hamiltonian 
is  taken  as  the  sum  of  Dirac  Hamiltonians,  plus  the  usual  Goulombic 
repulsion  terms  between  electrons.  This  Hamiltonian  is  thus  employed 
in  a  matrix  SCF  method  using  as  basis  functions  Slater-Type-Orbitals 
with  non-integer  principle  quantum  numbers.  The  usual  Breit  inter¬ 
action  operator  is  used  as  the  relativistic  interaction  term,  and  is 
treated  in  first-order  perturbation.  The  program  constructed  using 
this  method  was  employed  to  make  relativistic  calculations  on  He, 

Be,  and  Ne  (all  1Sq) . 

Several  calculations  considerably  beyond  Hartree-Fock  approxi¬ 
mation  have  also  been  carried  out.  Thus  "correlated"  orbital  calcu¬ 
lations  have  been  carried  out  for  He-like  systems  and,  in  a  most 
comprehensive  series  of  calculations,  the  ground-state  wave  functions 
of  He,  Li,  and  Be  were  approximated  by  a  superposition  of  configura¬ 
tions  with  35  terms  for  He  and  up  to  55  terms  for  Be.  Isoelectronic 
series  were  also  studied  for  this  series  of  atoms. 

^Theoretical  work  on  atoms  is  also  included  here  because  of  the  fre¬ 
quent  role  atomic  calculations  have  played  as  a  preliminary  to  the 
corresponding  molecular  calculations,  directly  or  indirectly. 
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B.  VERY  ACCURATE  CALCULATIONS  ON  THE  HYDROGEN  MOLECULE 

A  number  of  Important  results  for  the  ground  and  certain  excited 

states  of  Hg  have  been  obtained.  These  calculations  employed  a  James - 

Coolidge  type  expansion  to  represent  the  electronic  wave  function. 

Initial  calculations  for  H2(X  1Sg)  employed  up  to  a  40-term  expansion 

and  together  with  a  Hartree-Fock  wave  function  gave  the  R  variation 

of  the  correlation"  energy.  These  calculations  were  later  extended 

to  permit  up  to  80  terms  in  the  expansion,  and  the  James-Coolidge 

technique  was  generalized  to  permit  the  correct  asymptotic  behavior 

at  large  R.  Thus  calculations  were  made  on  the  X  12+.  b  c  1n 

1  +  g  uJ  u' 

and  B  2U  states  of  H2  and  the  first  two  12+  states  of  HeH+.  The 

resulting  wave  functions  were  then  used  to  obtain  nuclear  motion  and 
relativistic  corrections  and  the  final  energies  are  within  a  few 
kayers  (cm*1)  of  experiment  at  Re.  These  wave  functions  have  been 
extensively  used  to  calculate  a  number  of  expectation  values  and 
molecular  properties;  e.g.  the  parallel  and  perpendicular  polariza¬ 
bility  of  hydrogen,  the  dipole  moment  of  HD  etc.,  all  at  several- 
internuclear  separations.  In  a  parallel  study  of  Hg(X  1z+) ,  the 
complete  4-particle  problem  was  variationally  studied,  l.e;  without 
assuming  the  Born-Oppenheimer  separation.  A  number  of  related  studies 
were  also  completed  and  this  general  program  and  certain  ramifications 


continues . 
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C.  HARTREE-FOCK  AND  OTHER  CALCULATIONS  FOR  DIATOMIC  MOLECULES 

Calculations  cn  diatomic  molecules  have  usually  dominated  the 
computational  efforts  and  are  often  the  results  of  straightforward 
extension  to  diatomic  systems  of  theory  developed  for  atoms.  At 
the  beginning  of  this  period,  one  generation  of  diatomic  calcula¬ 
tions  was  Just  concluding;  during  this  period  the  second  generation 
of  calculations  was  organized,  programmed,  extensively  carried  out, 
and  now  at  the  conclusion  of  the  period,  a  new  "Generalized" 
Hartree-Fock  procedure  is  being  organized  and  programmed,  —  the 
third  generation. 

The  "first"  generation  of  calculations  were  minimal  basis  set 
LCAO-MO-SCF  calculations  for  a  number  of  diatomic  molecules  with 
or  without  optimization  of  orbital  exponents.  Thus  results  were 
obtained  for  Lig,  Beg,  Cg,  Ng,  Fg,  LiH,  BH,  NH,  HF,  CO,  BF,  and 
LiF,  all  for  states  (which  is  not  the  ground  state  for  Cg 

and  NH) .  For  these  molecules,  certain  population  analyses  were 
carried  out  and  with  a  calculated  potential  curve,  spectroscopic 
constants  were  obtained.  For  this  same  sequence  of  molecules, 
limited  configuration  interaction  calculations  were  also  carried 
out  using  the  minimal  basis  set  results  already  obtained.  These 
results  were  not  to  Hartree-Fock  accuracy  and  suggested  more  ex¬ 
tensive  calculations  were  needed. 

Between  these  results  and  the  "second"  generation  were  several 
double-zeta  LCAO-MO-SCF  calculations;  thus  slightly  more  extended 
calculations  were  carried  out  for  Ng  and  N+  and  for  Og.  Finally, 
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a  relatively  extended  basis  set  calculation  was  made  for  HF,  which 
was  the  closest  yet  to  a  molecular  Hartree-Fock  result  (except 
for  Hg) , 

The  "second”  generation  of  diatomic  calculations  began  with 
the  completion  of  two  completely  revised  and  extended  (e.g.  to  study 
open-sheli.  cases)  machine  programs;  one  for  homonuclear  and  one  for 
heteronuclear  diatomic  molecules.  This  also  coincided  with  a  new 
generation  of  high-speed  computing  equipment.  Two  to  three  years 
of  heavy  computation  followed  and  continues.  These  new  codes  and 
Increased  machine  speed  permitted  optimization  of  orbital  exponents 
and,  more  importantly,  allowed  large  enough  basis  sets  to  obtain 
true  Hartree-Fock  orbitals  for  many  diatomic  molecules  using  the 
matrix  method  of  Roothaan.  Thus  for  a  selected  group  of  molecules, 
convergence  studies  were  made  to  establish  the  necessary  basis  set 
size,  composition,  and  extent  of  optimization  of  orbital  exponents. 

This  ground  work  laid,  a  very  broad  and  Intensive  series  of  calcu¬ 
lations  were  started.  These  calculations  are  outlined  as  follows: 

1 .  Results  for  First  Row  Homonuclear  Diatomic  Molecules,  .  — 

Calculations  have  been  completed  (or  near  completion)  for  Lip(X  12+) , 
LI+(X  22g,  A  2nu),  Be2(V),  n2(X  5S~) ,  Cg(a  X2+,  Af  V),  Ng(X  , 
N+(X  2Z+,  A  B  2Z+),  02(X  h-g,  a  ^  b  h+g) ,  0+{X  Hg) ,  0*(X  2Hg) , 

F2(X  ^2g)  and  F^  Ions.  Most  of  these  results  Include  a  calculation 
of  potential  curves  and  a  detailed  calculation  of  molecular  property 
functions,  M(R) .  Gradual  extension  of  the  scope  of  these  investiga¬ 
tions  has  delayed  their  preparation  for  publication.  These  results 
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are  believed  to  be  quite  good  approximations  to  the  true  Hartree- 
Fock  orbitals  and  results.  Much  less  accurate,  minimal  LCAO-MO-SCF 
results  nave  been  obtained  for  the  second  row  homonuclear  series, 

Na0  to 

d  CL 

2.  Results  for  First  and  Second  Row  Diatomic  Hydrides,  AH.— 
Wavefunctlons  to  Hartree-Fock  accuracy  have  been  obtained  for 
L1H(X  1S'1'),  BeH(X  2Z+) ,  BH(X  1Z+) ,  CH(X  2nr) ,  I,H(X  5Z') ,  OH(X  2n±) , 
HF(X  12+),  Ne.H(X  12+) ,  MgH(X  2Z+) ,  AiH(X  12+) ,  S1H(X  2nr) ,  PH(X  \~) , 
SH(X  2I1^)  and  HCi(X  ^Z+)  at  Re(Exptl.),  and  as  a  function  of  Inter- 
nuclear  separation.  These  results  Include  a  study  of  convergence 
to  the  Hartree-Fock  limit.  From  the  potential  curves,  spectroscopic 
constants  were  calculated  and  compare  very  favorably  with  experiment. 
These  wavefunctlons  have  also  been  extensively  employed  to  calculate 
molecular  properties  as  a  function  of  R.  The  "correlation'’  energy 
of  the  hydrides  is  discussed  for  these  sequences  of  hydrides  in 
comparison  with  corresponding  separated  atoms  and  united  atom  values 
appropriate. 

These  ground  state  AH  results  have  also  lead  to  preliminary, 
but  very  good,  wavefunctlons  for  the  corresponding  AH+  series,  a 
number  of  excited  states,  and  a  few  negative  ions.  For  example, 
calculations  are  started  on  the  excited  states  of  BeH(lC22C2l7r,  A  2n) 
and  MgH(KL4c  27t,  A  n)  to  eventually  calculate  the  transition  proba¬ 
bility  for  the  X  2Z+— A  2n  transition.  Calculations  were  also  started 
for  the  0H(A  CZ+)  and  SH(A  ^Z+)  excited  states  also.  A  Hartree-Fock- 
Roothaan  wavefunctlon  (with  optimized  orbital  exponents)  was  also 
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obtained  for  OH'fX  V)  for  various  R  values  which  permitted  an 
estimate  of  the  spectroscop'  ;  constants  of  OH',  which  is  in  very 
good  agreement  with  experiment. 


3-  Miscellaneous  Results  for  Other  Hetermnn.^i,^  Diatom^ 
"oleoules.-  The  most  extensive  non-hydride  heteronuclear  results 
have  been  for  CO(X  V)  and  BF(X  V)  which  are  near  comparable 
accuracy  to  the  results  for  N2(X  .  A  large  number  of  calcula¬ 
tions  of  a  preliminary  nature  were  carried  out  to  obtain  a  reposi¬ 
tory  of  wavefunctions  for  many  heteronuclear  diatomic  molecules. 
These  calculations  made  use  of  Information  already  gained  in 
choice  of  basis  set  composition  and  orbital  exponents.  Results 
for  approximate  50  molecules  and  states  were  thus  obtained,  although 
in  most  cases  Hartree-Fock  accuracy  is  not  claimed.  These  include 
results  for  one  13  electron  sequence  PeF,  BF+,  BO,  C0+,  but  not  ON, 
each  for  the  X  2  ,  A  1^,  H  2nr  states,  results  for  BeO,  DO,  DP, 
LICf,  NO,  and  others.  Among  other  things,  these  calculations 
suggest  a  \-h  band  system  for  DO  in  the  red  or  Infrared  Experi¬ 
mentally  unknown)  and  an  electron  affinity  for  BO  of  about  1  eV. 

The  first  beginnings  of  the  third  generation  ol  diatomic  calcu¬ 
lations,  again  coinciding  with  a  new  generation  of  computers,  is 
the  development  and  demonstration  of  the  value  of  the  "Generalized 
Hartree-Pock"  method.  In  this  approach,  a  relatively  small  number 
of  configurations  are  mixed  with  the  ground  configuration,  and  the 
configuration  mixing"  problem  is  solved.  However,  the  new  orbitals 
which  give  rise  to  the  excited  configurations  are  also  optimally 
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determined  in  concert  with  the  ground  state  or  "core"  orbitals. 

Thus  a  doubly  iterative  method  obtains  as  the  configuration  mixing 
coefficients  and  all  orbitals  are  successively  re-determined  by  the 
variational  principle  to  a  certain  threshold  for  both  configuration 
mixing  coefficients  and  molecular  orbitals.  By  a  judicious  choice 
of  the  excited  configurations  as  those  arising  only  from  exciting 
"valency  shell"  molecular  orbitals,  the  calculated  potential  curves 
goto  the  correct  values  at  large  R  (in  contrast  to  the  Hartree-Fock 
curves)  and,  in  the  calculations  carried  out  to  date,  the  potential 
curves  calculated  are  in  very  good  agreement  with  experiment. 
Generalized  Hartree-Fock  calculations  were  carried  out  for  Hg,  Lig, 
and  Fg  with  greater  emphasis  on  Lig.  These  calculations  served  as 
an  important  trial  case  for  examining  the  feasibility  and  accuracy 
of  the  Generalized  Hartree-Fock  method  and  the  successes  have 
prompted  major  programming  efforts  to  extend  these  calculations  to 
diatomic  molecules  in  a  general  way. 

Certain  other  miscellaneous  diatomic  calculations  have  also 
been  carried  out.  These  include  the  extensive  configuration  in¬ 
teraction  results  for  the  he-He  interaction  potential  to  resolve 
an  experimental-theoretical  disagreement  and  a  very  illustrative 
configuration  Interaction  calculation  for  Hg(X  V)  to  study  the 

correlation  energy  as  associated  with  specific  kinds  of  configura¬ 
tions. 
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D.  CALCULATIONS  ON  POLYATOMIC  MOLECULES 

The  uaxealations  on  polyatomic  molecules  Involved  limited 
LCAO-MO-SCF  calculations  on  certain  closed  shell  linear  molecules 
and  the  one-center  expansion  SCF  calculations.  The  LCAO-MO-SCF 
calculations  were  an  outgrowth  of  the  earlier  minimal  basis  set 
calculations  on  diatomic  molecules,  and  were  not  near  Hartree- 
Fock  accuracy.  These  calculations  were  for  C02,  C2H2,  C7,  C^,  N-^, 
N02,  HF2,  HCN,  and  (CN)2;  all  linear  molecules,  each  at  the  equi¬ 
librium  configuration  (if  known) .  In  several  cas.s,  double-zeta 
basis  sets  were  employed  (HF2  and  CO^) .  These  calculations  were 
often  accompanied  by  results  for  certain  molecular  properties  and 
interpretative  data,  and  have  served  as  a  useful  starting  point 
for  further  work  which  is  in  course  elsewhere.  In  another  minimal 
LCAO-MO-SCF  calculation,  CH^  was  studied  including  all  integrals. 

The  one-center  expansion  SCF  calculations  involve  no  basis 

functions  on  the  peripheral  atoms,  always  hydrogen,  for  molecules 

AH  .  The  initial  set  of  one -center  calculations  did  not  use  a 
n 

symmetry-specific  program,  which  limited  the  basis  set  size,  and 
calculations  were  carried  out  for  HF,  CH^,  SiH^,  NH^,  NH^,  PH^,  Ph£, 
HgO,  HgS,  and  HCi .  These  calculations  clearly  indicated  for  the 
first  time  the  limitations  and  possibilities  of  the  one-center 
expansion  SCF  method.  Crude  but  useful  wavefunctions  were  obtained 
for  certain  modera"1  ely  large  polyatomic  molecules.  The  second 
phase  of  the  one-center  SCF  calculations  involved  a  s’  metry- specific 
program  for  C^v  and  C2v  symmetry.  Using  these  programs,  extensive 
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one-center  SCP  calculations  were  carried  out  for  NH^  to  study 
the  inversion  barrier.  Calculations  for  the  isoelec tronic  mole¬ 
cules  OH^  and  CHg  have  also  been  completed.  The  eight  electron 
sequence  BeH^,  BH^,  CHl,  NH*  and  OH were  also  considered  and  all 
are  predicted  to  be  planar  in  their  ground  states.  In  all  of 
these  later  calculations,  very  large  basis  sets  (including  d,  f, 
g,  and  h  STO’s)  have  been  employed  and  the  orbital  exponents  were 
carefully  optimized.  A  similar  computer  program  was  c  mpleted  for 
C2v  symmetry,  and  some  preliminary  calculations  were  mad?  for  HgO. 
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